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Abstract 

We discuss the pure gauge Schwinger-Dyson equation for the gluon propagator in 
the Landau gauge within an approximation proposed by Mandelstam many years 
ago. We show that a dynamical gluon mass arises as a solution. This solution is 
obtained numerically in the full range of momenta that we have considered without 
the introduction of any ansatz or asymptotic expression in the infrared region. 
The vertex function that we use follows a prescription formulated by Cornwall to 
determine the existence of a dynamical gluon mass in the light cone gauge. The 
renormalization procedure differs from the one proposed by Mandelstam and allows 
for the possibility of a dynamical gluon mass. Some of the properties of this solution, 
such as its dependence on Aqcd and its perturbative scaling behavior are also 
discussed. 
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1 Introduction 



It is widely believed that Quantum Chromodynamics (QCD) is the theory 
which describes the strong interaction. For this theory we know that pertur- 
bation theory has become a reliable field theoretical method of calculating 
and predicting most of the quantities in processes where high energies are 
transferred between quarks and gluons. 
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This successful procedure to deal with the strongly interacting phenomena is 
known to be inadequate when it is applied to the infrared region. There are 
phenomena at low energies such as dynamical chiral symmetry breaking, that, 
in principle, could only be described when all orders of perturbation theory are 
taken into account, it means that they are necessarily of a non-perturbative 
nature. 

To bridge the gap between these two regions, infrared and ultraviolet, two main 
non-perturbative approaches are available, the lattice theory which is based 
on discretization of space-time and a continuum one which makes use of an 
infinite tower of coupled integral equations that contain all the information 
about the theory - the so called Schwinger- Dyson Equations (SDE). 

In the continuum, we hope that the SDE provide an appropriate framework to 
study the transition from the perturbative to the non-perturbative behavior 
of the QCD Green functions. However, its intricate structure only become 
tractable when we make some approximations. 

Many attempts have been made to understand the gluon propagator behavior 
through SDE. In the late seventies Mandelstam initiated the study of the 
gluon SDE in the Landau gauge [1]. Neglecting the ghost fields contribution 
and imposing cancellations of certain terms in the gluon polarization tensor, 
he found a highly singular gluon propagator in the infrared. This enhanced 
gluon propagator was appraised for many years in the literature, firstly because 
it provided a simple picture of quark confinement [2], since it is possible to 
derive from it an interquark potential that rises linearly with the separation, 
and secondly because a gluon propagator, which is singular as 1/ g^, has enough 
strength to support dynamical chiral symmetry breaking, as it was claimed 
in the studies of the quark-SDE. This approximation and its solution were 
extensively studied by Pennington and collaborators [3]. However, these results 
are discarded by simulations of QCD on the lattice at 95% confidence level 
[4], where it is shown that the gluon propagator is probably infrared finite. 

Infrared finite solutions are also found in the Schwinger-Dyson approach, as 
result of different procedures. Many years ago, making use of the "pinch tech- 
nique", Cornwall built up a gluon equation trading the conventional gauge- 
dependent SDE by one formed by gauge-independent blocks. Analyzing this 
new equation, he obtained a gluon propagator endowed with a dynamical mass 
[5]. 

Recently, a quite extensive work on pure gauge SDE has been done by the 
authors of Ref. [6] where they have shown that when the ghost fields are taken 
into account, the gluon propagator is suppressed and the ghost propagator 
is enhanced in the infrared region. Such solution was shown to satisfy the 
Kugo-Ojima confinement criterion [7,8]. These propagators exhibit an infrared 
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asymptotic power law behavior which is characterized by a critical exponent k. 
Axiomatic considerations [9] and the latest results of the coupled gluon-ghost 
SDE seems to suggest that k = 0.5 is allowed [10], signaling the possibility of 
dynamical mass generation for the gluon. 



All these solutions appear because different approximations were used and 
furthermore it is also perfectly possible that in the same approximation more 
than one solution arise. It is interesting to note that, according to Mandel- 
stam's work (see the comments after Eq.(2.16) of Ref.[l]), a massive gluon 
solution was discarded from the beginning in his study. 



It is important to stress at this point that a dynamical gluon mass does not 
break gauge invariance and is consistent with massive Slavnov- Taylor identities 
[5] . We would like to emphasize that the presence of a dynamically generated 
mass also does not mean that gluons can be considered as massive asymptotic 
states similar to dynamically generated quark mass does not mean that quarks 
can be observed as massive asymptotic states. Why quarks and gluons are 
not observed as free states is the well known problem of confinement. In the 
case of a theory with massive gluons we know that such theories admit a 
vortex solution that may give a clue about the confinement mechanism [5,11]. 
Furthermore, a dynamically generated gluon mass is possibly connected to 
the existence of a QCD infrared fixed point [12], whose presence has many 
phenomenological implications as nicely reviewed in Ref.[13]. 



Our aim here is to revisit the gluon SDE within the Mandelstam approxi- 
mation, in order to obtain a massive gluon solution. In section II, we start 
building up the gluon SDE, which embodies not just the full gluon propagator 
but also involves the full triple gluon vertex. In order to allow that a dynamical 
gluon mass takes place without breaking the relationship between the Green's 
functions of different orders which are imposed by the gauge invariance, the 
full triple gluon vertex behavior is modeled by a suitable Slavnov- Taylor iden- 
tity. Having found the gluon-SDE, the ultraviolet divergences will be removed 
by a subtractive renormalization procedure, in a similar way as performed by 
Cornwall in Ref. [5] (although the equation is different), which introduces an 
arbitrary scale /x^ that can be related to the usual QCD scale Aq^d- This 
discussion will appear in section III. We then solve the gluon equation by an 
iterative numerical procedure on the whole range of momenta and present our 
results in section IV. In the present work we do not need to make any ansatz 
for the infrared solution. Section V contains a discussion about the vacuum 
energy and stability of the solution. We draw our conclusions in section VI. 
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2 The gluon equation in the Mandelstam's approximation 



The SDE are coupled integral equations which relate all the Green's functions 
of the theory. To illustrate how intricate is its structure we can look at what are 
the Green functions which are involved in the full gluon equation. Neglecting 
the fermionic interactions we can see, in Fig.(l), that the gluon propagator is 
written in terms of itself, the full 3 and 4-point gluon vertex, F^^^p and T ^^p^, 
and also the full ghost propagator and the gluon-ghost coupling. 




Fig. 1. The complete Schwinger-Dyson equation for the gluon propagator without 
quarks. The spiral lines represent gluon field and the dashed lines ghosts. The black 
blobs are the full propagators while the white ones are the full vertices. 

Actually, these unknown three and four point-functions obey their own SDE, 
which involve higher n-point functions which naturally must satisfy, in their 
turn, their own SDE. In fact, it is this entanglement of equations which makes 
unavoidable the use of some truncation schemes. 

One famous truncation scheme is the Mandelstam's approximation where the 
fermion fields are neglected, since we believe that a pure Yang-Mills theory 
must carry all the main features of QCD. Furthermore, based on perturbative 
results, the ghost fields are also neglected. The justification, for the latter 
approximation is that the contribution which comes from the ghost field is 
supposed to be small even in the non-perturbative region. Therefore, this 
approximation can be represented pictorially by the graphics which compound 
the first line of the Fig.(l) and is written as [1] 
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where p — k + q, Dq and Fq are respectively the propagator and three- gluon 
vertex at tree level, while D and F are the two and three points full Green 
functions. In the case of the full gluon propagator in Landau gauge, D^"^, we 
can write it as 

where it will be useful to define the function, D{q^), in terms of the gluon 
renormalization function, Z[q^). 

Diq^) = (3) 



Neglecting all contributions coming from ghosts fields, the Slavnov-Taylor 
identity between the three gluon vertex and the inverse of the gluon prop- 
agator can be expressed in the following simple form 



In order to allow that a massive gluon propagator will be also compatible 
with the Slavnov-Taylor identity expressed above, and supposing that the 
gluon renormalization function, Z[q^) admits an expression of the following 
form 

^(?') = (5) 



we note that it must be added new terms, that have massless poles, to the 
structure of the three gluon vertex, which, apart from a group theoretical fac- 
tor, lead us to a first modification that we should introduce in the construction 
of the full vertex, which is the one prescribed by Cornwall many years ago [5] 



r&(A:,P, q)-{k- p),5,. + ^^^^^^^1^ + c.p., (6) 



where c.p. means cyclic permutation, and, as discussed in Ref . [5] , Tq^^^ is the 
vertex for the massive theory. 
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As remarked by Corwnall [5], during the procedure of construction of a ver- 
tex function which automatically satisfied the Slavnov- Taylor identities, Ball 
and Chiu [14] have done a crucial assumption in order to get a unique form 
for the longitudinal vertex. They supposed that the vertex should be free of 
kinematic singularities, since most of these singularities violate the general 
analyticity requirements of the vertex. Therefore, tensorial structures which 
have massless poles in three gluon vertex were explicitly excluded in their 
construction, despite the fact that they have already mentioned in the same 
work, that there are certain types which might naturally occur without the 
breakdown of analyticity in the three gluon vertex, and one structure of this 
type is the one given by the second term in the right hand side of Eq.(6). 

It is important to bear in mind that the mass m in Eq. (6) has a momentum de- 
pendence and do not destroy the unitary behavior of the theory [5] . Moreover, 
as stressed before, the role of the latter term from Eq.(6), is only to allow that 
the gluon propagator could assume a non-zero value, i.e. D^^{q^ = 0) 7^ 0, 
at the deep infrared region, without breaking the gauge invariance imposed 
by the Slavnov- Taylor identity. Exactly the massless poles of Eq.(6) lead to 
the possibility of a mass gap, contrarily to the solutions for the infrared gluon 
propagator behaving as found in Ref. [1,15] where a different vertex choice 
is made. 



The next step in the Mandelstam approximation is to define the final expres- 
sion for the full three gluon vertex. The form of the full three gluon vertex is 

[3] 

T^'^'ik.P.q) = -^^Tlr^'"''{k,p,q). (7) 



The use of the above full vertex, which is a combination of the Cornwall's and 
Mandelstam's prescriptions, simplify even more the structure of gluon SDE 
than the use of the bare vertex, once it implies a cancellation between the 
gluon renormalization function which comes from a full gluon propagator and 
the one that comes from the full triple vertex. 




Fig. 2. The gluon Schwinger-Dyson equation in the Mandelstam's approximation. 
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With the above expressions for the gluon propagators and vertices and con- 
tracting this result with the projector, proposed by Brown and Pennington 
[3], 

n^-{k) = 5^"" - (8) 



the follow equation comes out, 



gl 'Jdq^ [7q' 17 9\ 








167r2J 

fc2 




2{q% (9) 



where go is the bare coupling, we use the color factor, C2 = 3 and A is an 

ultraviolet cutoff which was introduced in order to render the integral finite. 
In fact, the great advantage of the Mandelstam approximation is that all the 
angular integrals can be performed analytically leading us to a much simpler 
equation. 

A remark concerning the massive term of the three gluon vertex is in order. 
This equation is exactly the same one obtained in Ref.[3]. Although the three 
gluon vertex has an extra term, no contributions come from it. Therefore, as 
mentioned before, its only role is to allow that the inverse gluon propagator 
can be different from zero at the origin consistently with the Slavnov-Taylor 
identity. 



3 Renormalization 



The procedure to perform the renormalization of SDE in the Mandelstam's 
approximation consists in introducing the gluon and vertex renormalization 
constants, Z3 and Zg, respectively which will absorb ultraviolet divergences of 
the equation. Through these constants we can define the following renormal- 
ized quantities 



D(q',A')^Zsi^i',A')DR(q',^i') 

go(A')^Zg(^i',A')g(^i% (10) 

where Dji and g are the renormalized gluon propagator and the renormalized 
coupling. 
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Substituting into the SDE, Eq.(9), the nonrenormahzed quantities, D and qq, 
by the renormahzed ones lead us to 

(11) 

where is given by 




The vertex renormahzation constant Zg can be ehminated from the Eq.(ll), 
using the identity, ZgZ^ — 1, which is only valid in this approximation [16]. 
However, there still remains the gluon renormahzation constant Z^ that is a 
naturally divergent quantity. 

Despite all efforts to obtain a totally consistent renormahzation of the gluonic 
SDE, be it in this truncation or even beyond the Mandelstam approximation, 
we know that the renormahzation of SDE is highly nontrivial. To illustrate this 
claim we notice that only in the nineties, Curtis and Pennington pointed out 
a truncation scheme, for QED, which is gauge independent and also respect 
the multiplicative renormalizability [17]. In the following we consider Z^ as 
a factor which renders the divergent terms of the gluon SDE in a finite one 
through a subtr active renormahzation. 

After imposing ZgZ^ = 1 in Eq.(ll) we are left with the following equation 



ez^ + zizi 



9^ 
167r2 



Z3 + 



A2 



- / dq' (-7) Dn(q') 

fe2 



+ 



+ 



47r 
47r 



'7q^ 17 q^ 9^ 





A2 



fc2 



(13) 



where q;(//^) = ^ 'f ^ ■ Z^ behaves as 1 plus an infinite piece. The factor 1 
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is consistent with the pcrturbative behavior, and the infinite part of Z3 must 
cancel the infinite part that comes out fi^om the first integral in Eq.(13), which 
is the only divergent integral in the above equation. 

It is important to notice that Eq.(13) is not of the same form as the one shown 
in Ref. [5] , as well as we shall not follow the same procedure adopted up to now 
in the many papers to renormalize the Mandelstam equation. In the previous 
papers where this equation was solved it was verified that the dressing of the 
gluon propagator could behave as 



The numerical solutions were then obtained assuming A — and B — 1, which 
are not independent conditions. Note that these conditions are necessary to 
obtain a solution compatible with the Slavnov- Taylor (ST) identities in the 
case of a massless gluon. As discussed in the previous section, it is possible 
to introduce a new piece in the gluon vertex that makes the ST identities 
compatible with dynamical mass generation for the gluons, alleviating any 
condition on A and B. 

We will face Z3 as the term that will eliminate the infinite contribution of the 
following expression 

Zs-^l dq' 7Dn{q') . (15) 



Of course, in the above subtraction there remains a finite contribution. Dii{q^) 
must be known to perform the above calculation, but it is going to be known 
only after we impose the renormalization conditions 

Mf^') = 1, (16) 



or equivalently 

Dn\n')=n'. (17) 



The role of the above equation is to guarantee that exactly at the point jj,'^ 
we recover the bare pcrturbative behavior, DR^q^) = 1/q^. For this reason it 
is important that the value of /i^ be fixed at a typical pcrturbative scale, i. e. 
it must satisfy the condition /i^ ^ ^qcd- 

Let us now discuss the prescription for the finite contribution to Eq.(15). The 
solution usually found for the Mandelstam gluon propagator is of the form 
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Di{k!^) ^ 1/k'^. But in principle, as long as wc do not impose A = and 
i? = 1 as discussed after Eq.(14), we cannot discard a solution behaving like 
D2{k^) ~ + ln(A;^ + m^)], which is a massive solution similar to the 

one found by Cornwall in Ref.[5]. 

If DR{q^) behaves as Di{q^) it is natural to have a finite contribution like 
Ki/k'^ in the integration of Eq.(15). The same is true for high values of k"^ 
if D^i^q'^) behaves as 1)2(9^), since the integral that appears in Eq.(15) will 
generate terms like InlnA^, Inln/c^ and again a term proportional to m^/A;^. 
It is also not hkely that the integral in Eq.(15) will give finite contributions 
hke 5n/{k'^Y n > 2, because this behavior is not consistent with any 
expected form of DR^q^). As ki and have exactly the same dimension, the 
simplest form of Eq.(15)which is compatible with both behaviors is 



where n has squared mass dimension and is going to be fixed by the renor- 
mahzation procedure = /i^). 

A more complex expression on the right-hand side would introduce new con- 
stants besides k, which could be fixed only in an ad hoc fashion. The 1 in the 
right-hand side of the above expression corresponds to the first Z3 term. The 
choice we made in Eq.(18) does not obey the conditions A — and B — 1 
that we discussed before, and, in principle, it allows even for a massive gluon 
solution. This possibility is an actual one as long as we obtain a stable solution 
for the final equation. 

In all the procedures used up to now to solve the Mandelstam equation, it 
has been assumed a cancellation of certain terms, and there is not any discus- 
sion about the finite terms that survive renormalization. It is known that in 
perturbativc calculations such finite terms are nothing else than small quan- 
tum corrections which modifies the tree level value of the renormalizcd quan- 
tity. On the other hand the physical quantities which are related to the non- 
perturbative dynamical mass generation will be exclusively generated by the 
quantum corrections. This means that such finite terms surviving renormal- 
ization are fundamental to the final result. 

With the definition described by the Eq.(18) and imposing the renormalization 
condition (17) into Eq.(13) we can obtain the following expression for the 
parameter 




(18) 
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and, consequently, the final expression for the inverse of the gluon propagator 
can be written as 




Since we are now dealing only with renormalized quantities, in the sequence 
we will dismiss the subscript R in the Green functions in order to get a more 
compact notation. It is interesting to note that the renormalization constant 
as it is established by Eq.(18) is proportional to 1 plus a function of /^^ and 
A^. Such behavior is compatible with the expected weak coupling expansion 
for this constant. 

As already stated, the renormalization procedure in Schwinger-Dyson equa- 
tions has an intricate structure and other choices have already been applied 
in these equations, such as the "plus prescription" [18] where the contribu- 
tions which violate the massless Slavnov- Taylor identity would be subtracted 
out of the right hand side of Eq.(9) or even in more elaborated CclSGS, clS the 
gluon-ghost coupled system, where we have to deal with a bigger number of 
renormalization constants, different approaches can be considered [6,19]. As 
long as we do not have an exact procedure for the renormalization problem 
in non-Abelian SDE, we have to face this prescription as one more try, that 
certainly can be improved, where the quantitative perturbative behavior of 
the renormalization constant Z3 is reproduced. 

Our prescription in Eq.(18) acts as a seed for the dynamical gluon mass gen- 
eration, since k defined by Eq.(19) is precisely a constant, like the A term in 
Eq.(14). Of course, this is not the most general prescription, and it does not 
reproduce the exact perturbative ultraviolet behavior of the gluon propagator. 
This behavior can be obtained if we use the following prescription 
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It is clear that the last term of Eq.(21) arises naturally from the integration of 
the perturbative propagator {D{q^) = in Eq.(15). However, the new con- 

stant (^2) cannot be fixed, since we have only one renormalization condition. 
In the next section we shall discuss the effect of a term like this one. 

Before starting the numerical calculation of the gluon propagator, it is use- 
ful to remind the definition of the running coupling that will be valid from 
the non-perturbative to the perturbative region. Remembering that, in the 
Mandelstam approximation, Z2,Zg — 1, it follows from Eq.(lO) that 



where we assume that Z{ii^) = 1. 

It is important to keep in mind that, in QCD, we have the possibility to define 
the running coupling constant using the different vertices of the theory, such 
as, the three gluon, four gluon, ghost-gluon or even the fermion-gluon vertices. 
Despite the fact that in the perturbative QCD all these constructions, based on 
different vertices, converge to the same behavior for the running coupling, since 
the Slavnov-Taylor identities are preserved, the same does not happen in the 
low energy scale, due to the complex renormalization procedure in the infrared 
region of QCD, which cause the loss of the multiplicative renormalizability. As 
a consequence we may not have a unique definition for the running coupling 
in the infrared region. For this reason it is important to stress that, in this 
approximation, our QCD running coupling is constructed on the basis of three- 
gluon vertex, since all the others vertices do not appear in this approximation. 

In the sequence we shall need the perturbative expression for the running 
coupling, which is given by 



where Aqqd is the usual scale for QCD and (3o is the first coefficient of the 
the Callan-Symanzik (3 function, (3{g) = ii{dg/dii), which is written as 



(22) 




(23) 




-A 



(167r2) 




+ ..., 



(24) 



where the (5q and /?i coefficients are given by 



/3o = ll - g% 
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/3i = 102--n;. 



(25) 



Since we neglected the fermions fields, we have that n/ = which lead us to 
/5o = 11 at one loop. 



4 The mcissive solution 

It has been shown in several works that it is not an easy task to obtain an 
analytical solution for the Eq.(20) [1,3,6] and therefore the only possibility to 
face this problem is through a numerical approach. 

Our aim here is to find consistent solutions over the whole momentum range 
without impose, neither in the infrared region nor in the ultraviolet, any pre- 
vious asymptotic behavior obtained from an expansion of the gluon rcnormal- 
ization function, Z[q^), at small (f or obtained from the known perturbative 
behavior. 

For this reason, we apply for the integral equation, given by Eq.(20), an itera- 
tive numerical method, starting with a trial function which can be too remote 
from the exact solution. 

To implement so, it is convenient to introduce the following variables x — k'^ 
and y = in the Eq.(20), and we also recall that since the beginning of this 
article all variables are in the Euclidean space, with these changes we can 
write Eq.(20) as 

= « + . + A / * - ^ - yDiy) + a/ * yD(ym 

where 

^ - -V* (8-0-9 ^^(*) - (^) (^^' 

with A = Q;(/x^)/47r. 

In order to study, in more details, the small x region, wc use a logarithmic 
grid for the variables x and y, which allow us to vary the momentum from the 
deep infrared to the ultraviolet region. Such logarithmic grid split the whole 
momenta range in two regions: the infrared region - defined by the range 
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[0, //^] and the ultraviolet that comprehends the range [/i^,A^]. The aim of 
this separation is to allow us to set 2{ii^) = 1 or equivalently, D~^[^'^) — /s^ 
[21]. 

We start with a trial function D{x) and use a cubic spline interpolation for 
generating the values of D{y) which will be utilized in the right hand side of 
the Eq. (27) and in the sequence wc compute this integral through the Adaptive 
Richardson- Romberg extrapolation. 

The initial trial function is compared with the result which was obtained after 
the integration and the convergence criteria to stop the numerical code is to 
impose that the difference between input and output functions must be smaller 
than 10-^. 

It is important to stress that we have tested variations of initial guesses and 
verified that our results are completely independent of the trial function im- 
posed for D{x). 

We also analyzed the solution proposed in Ref.[3] and we noticed that it can 
only be reproduced if we consider exactly the same momentum range showed 
in their Fig. 3 (as is already stressed in Ref.[16]) and use a trial function which 
is very close to the result found in [3] . Such exigency reflects the instability of 
this solution, because if we extend the numerical range or even start from a 
different guess we do not recover the divergent behavior. 

Our input data arc the rcnormalization point, jj,^, and the coupling constant 
defined at this point, a{fi^), which has the effect of fixing the value of the 
QCD scale, Aqcd, through the Eq.(23). 

As stressed in the Sec. (3), the value of /i^ must be fixed at a typical pcrturba- 
tive scale in order to recover the high energy behavior of the gluon propagator. 
With the aim of analyze the dependency of our solution on the renormalization 
point, we vary the values of /x^ and a{fi^) within the range [10 GeV^, 30 GeV^] 
and [0.20, 0.25] respectively. Such variation correspond to run the Aqcd pa- 
rameter from 182 Me V to 557 Me V as we show in the Table(l). 

The curves produced by these different scales are plotted in the Fig. (3) where 
it is shown the gluon propagator, D[q^), as a function of the momentum 
q^. The external curves delimit the lower and the higher values of the Aqcd 
in the range mentioned above. The other set of input, shown in the Table(l), 
reproduce the same qualitative behavior and they are restricted to the shadow 
band. 

We have also computed the case where a{fi^) is fixed at the bottom quark 
mass, ml = (4.5)^ GeV^, and its central value is a{ml) = 0.22 [22], such 
solution is represented by the curve "fine -|- circle" displayed on the Fig. (3). 
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Fig. 3. Gluon propagator, -D(g^), as function of momentum for different scales. 

The line + square curve was obtained when a(10 GeV^) = 0.2 which corresponds 
^QCD = 182 MeV, while in the line + triangle curve, a(30 GeV^) = 0.25, which leads 
to Aqcd = 557 MeV. The shadowed area delimits the curves with Kqcd varying 
within the range [182 MeV, 557 MeV]. The central curve (line + circle) was obtained 
when we fix the renormalization point, /.t^, at bottom quark mass, = (4.5)^ GeV^ 
with the central value of C({ml) = 0.22. 

We have run our numerical code efficiently within a momenta range of twelve 
orders of magnitude where the typical momenta values which were utilized 
vary from 10~^ to 10^ GeV^. We set the renormalization point at ml, which is 
located approximately in the middle of our typical momenta range, and where 
the physical quantities can be certainly described by the perturbative theory. 

It is interesting to provide an analytic expression for the gluon propagator, in 
order to analyze the gluon mass values obtained in infrared region, when we 
set the renormalization point, /.i^, at different values. With this aim we fit our 
numerical data by an Euclidean massive propagator expressed by [5] 



An 



(28) 




where 




(29) 
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We have also considered the simpler expression 

D{q') = 2^ (30) 



where the dynamical mass Ai^{q^) is described by 

M\f) = (31) 



The last fit can be motivated by the gluon polarization tensor behavior at 
high energies, which can be predicted by OPE as [23] 



IIope{P ) ~ p-2 ■ (32) 



On the other hand recently there has been a lot of discussion about a possible 
bilinear condensate of the gluon field [24]. Such condensate, (A'^A^), would 
be responsible for a mass term appearing in the infrared gluon polarization 
tensor [24]. Therefore the fit provided by Eq.(31) is just the simplest way to 
account for the different condensate contributions to the polarization tensor, 
from where we could expect approximate relations of the form m'^ oc (A^A^) 
or m'^ oc i^^C^" G ^y"^ . We will discuss such type of relation in the next section. 

We can apply these simple fits for all curves shown in the Fig. (3) and in all 
cases it is remarkable the agreement found from the deep infrared region up 
to the ultraviolet regime, using a fit which has a unique parameter, mo or m' . 

In particular, we plot in the Fig. (4) the numerical solution for the gluon prop- 
agator. D{q'^), when /x^ = together with the curve obtained through our fit 
given by Eqs.(30) and (31) when m'^ = 1.11 GcV^ as well as the fit provided 
by the Eqs.(28) and (29) when mo^ = 0.82 GeVl 

In particular we can easily extract the value of m' for the others curves 
displayed in the Fig. (3), to do so we basically have just to note that in 
the limit of q^ —> 0, Eqs.(30) and (31) of our gluon propagator reduces to 
D(q^ — >■ 0) = 1/m'^, and therefore the inverse of m'^ is given by the value of 
the point in what the gluon propagator curves cross the y axis in Fig. (3). 

We verified that the value of itlq and m' depend on the choice of the renormal- 
ization point, /x^, however we must remember that when we change its value, 
actually what we are really changing is the scale of the theory, once that /^.^ 
and Aqcd are linked by Eq.(23). For this reason what matters is the analysis 
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Fig. 4. Numerical solution for the gluon propagator, D{q^), versus momentum 

for a{m,f^) = 0.22. Wc compare this numerical solution with the fit given by Eq.(30), 
where m'^ = 1.11 GeV^, as well as the one given by the Eqs.(28) and (29) when 
mo^ = 0.82 GeV^ 

of the ratio itiq/Aqcd ot m' /Aqcd which, in principle, give to us a better idea 
about the true dependency on the renormalization point of our solution. 

Fixing the coupling constant and running the renormalization point we can 
see from the Table(l) that the aforementioned ratios practically do not vary 
for the set of coupling constants shown in the table. This means that the ratios 
mo/ Aqcd ^'/Aqcd are quite stable in our procedure. 

Table 1 

Values of the renormalization point, /x^, and coupling constant, a(/Li^), used as input 
data in the Eq.(27). In the third column, we have the values of Aqcd computed 
with the usual perturbative value of /3o = 11. The values of the ratios mo /Aqcd 
and m' /Aqcd sue also shown in the last two columns. 



a(M^) Aqcd mo /Aqcd m' /Aqcd 

(/3 = 11) CornwaU OPE 



0.20 


lOGeV^ 


182 MeV 


3.18 


3.90 


0.20 


20GeV2 


257 MeV 


3.18 


3.90 


0.20 


SOGeV^ 


315 MeV 


3.18 


3.90 


0.22 


ml = 20.25 GeV^ 


335 MeV 


2.70 


3.15 


0.25 


lOGeV^ 


321 MeV 


2.23 


2.45 


0.25 


20GeV2 


455 MeV 


2.28 


2.44 


0.25 


SOGeV^ 


557 MeV 


2.23 


2.43 
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We have also checked the effect of the prescription shown in Eq.(21) attribut- 
ing random values to K2. We also obtained massive solutions with a better 
ultraviolet behavior. The numerical curves are totally identical to the ones in 
Fig. (3) but they scale up or down by a constant factor as we change the value 
of k2. The only way to eliminate such constant is forcing point by point of 
the numerical solution in the large momentum region to match with the curve 
given by the known asymptotic perturbative solution. 

We now turn to the infrared behavior of the coupling constant which was 
built based on the triple gluon vertex as outlined above. By imposing the 
limit ^ in the Eq.(22) where the rcnormalization function, Z[(^\ can be 
extracted from the fit which is expressed by Eq.(30) and (31) we clearly ob- 
tain a vanishing coupling in the deep infrared region within the Mandelstam 
approximation. Although it may look surprising to find that in a confining 
theory the IR coupling constant goes to zero, it is important to stress again 
that, for QCD, we have distinct definitions for which are based on dif- 

ferent vertices of the theory, moreover it is interesting to say that using the 
same vertex, a lattice QCD simulation found the same vanishing behavior for 
the couphng [25]. Of course, the introduction of ghosts can modify this result, 
and it seems unlikely then that this limit does refiect the true behavior of the 
coupling constant since we know that it must develop a non-trivial fixed point 
in QCD infrared region [12]. 

Following the same procedure that are perfomed here, the qualitative behavior 
of the gluon propagator does not change when the ghosts fields are included 
when we study the coupled SDE for the gluon and ghost [26]. We believe 
that the main role of the ghost fields in covariant gauges is to guarantee that 
the coupling constant will be finite and different from zero in the infrared 
regime, although it is quite hard to obtain numerically this freezing of the 
coupling without imposing any previous asymptotic form for the gluon and 
ghost propagators . 



5 Vacuum energy and stability of the solution 

In this section we would like to discuss some points about the stability of the 
massive solution in this approximation. We call attention to the fact that our 
solution is obtained numerically in the full range of momenta. Other solu- 
tions for this kind of equations in general assume one particular form of the 
solution in the far infrared region and adjust free parameters in a "in" and 
"out" procedure. We notice that our numerical code finds stability for the 
1/g*^ solution found in Ref. [3] only when we enter in the code a seed basically 
given by the final result of Ref. [3], otherwise there is no convergence up to 
a large computing time. We credit this behavior to the fact the l/gr^ must 
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be an unstable solution of the SDE in this approximation, because any seed 
that is shghtly away from the result of Ref. [3] does not lead to a convergent 
calculation. Unfortunately we cannot say more than that about the stability 
of the solution. 

It is possible to use some methods of integral equations to study the existence 
and stability of the solutions as in Ref. [27], but these methods may depend 
on the many approximations that are necessary to make in order to obtain 
a tractable equation before they can be apphed. We beheve that the best 
to be done is to compute the vacuum energy for composite operators [28]. 
This vacuum energy is a function of the full propagator and vertices of the 
theory [28,29]. The idea is simply to recall that the vacuum energy will select 
the solution that leads to the deepest minimum of energy as discussed in 
Ref. [30] in the case of pure gauge QCD. If we follow the results of Ref. [30] 
we can foresee that the massive solution is the one selected by the vacuum. 
However we can do more than that and we will show that the computation of 
the vacuum energy obtained in the previous section also leads to a consistent 
value for the gluon condensate. 

We will briefly outline the calculation of the vacuum energy with the for- 
malism of the effective potential for composite operators [28]. It will also be 
computed according to the Mandelstam's approximation, what is equivalent 
to neglect diagrams with fermions, ghosts and the quadrilinear gauge cou- 
phng. The details can be obtained in Ref. [31]. The effective potential with the 
approximations already discussed has the form [28] 

^^^^ ^ I / m'^^'^^'"^''^ ~ ^''^ + + ^^^^ 



where D{Dq) is the complete (bare) gluon propagator, ¥2(0) is a two-particle 
irreducible vacuum diagram and the equation 



gives the SDE for gauge bosons in the Mandelstam's approximation. 
The two-loop contribution to ¥2(0) is given by 

V2{D) = — rr(r(3)L»r(3)L»L»), (35) 



where T^^^ is the trilinear gauge boson coupling [30], and in Eq.(35) we have 
not written the gauge and Lorentz indices, as well as the momentum integrals. 
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The vacuum energy density is given by the effective potential calcnlated at 
minimum subtracted by its perturbative part, which does not contribute to 
dynamical mass generation [28,29] 

{n)^Vmin{D)-Vmin{Dp), (36) 



where Dp is the perturbative counterpart oi D. It is easy to see that [31] 



3{N^-l) f d'^P 



J (2^) 



n 



ln( 1 



n2 



p2 + n V PV 3P2(p2 + n) 



,(37) 



where all the quantities are in Euclidean space, N — 3 ior QCD and U is 
the gluon polarization tensor. We will compute Eq.(37) with 11 = A^^(P^), 
where AI^(P^) is given by Eq.(31) with m'^ = 0.50 GeV^, with the coupling 
constant, a{n'^ = lOGeV^) = 0.20. We obtain 

TT^ = 3.60. (38) 



with m'/AqcD = 3.90. 

Let us recall that the vacuum expectation value of the trace of the energy 
momentum tensor of QCD is [32] 

= ^ {G,.Gn , (39) 



where the perturbative P{g) function up to two loops is given by Eq.(24), 
where the coefficients Pq and Pi are expressed in Eq.(25), and with Q;(/i^) = 

We can relate the Eq.(39) to the vacuum energy through 

{^) = I (©mm) ■ (40) 



This last expression for (^2) can be compared with the value of Eq.(38) and 
in this way we obtain one estimative of the gluon condensate. 

Using Eq.(24) and Eq.(25) with Uf = 5 (assuming that the inclusion of 
fermions does not change drastically our results) and using a{fi^ = 10 GeV^) = 
0.20 (equivalent to m' / Aqcd = 3.90) we have found the following value 

—G^'^G^J) = 0.015 GeV^ (41) 
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We obtained this last value with Aqcd = 182 MeV in the left hand side of 
Eq.(38) consistent with the perturbative value shown in the third column of 
Table 1. The outcome of this procedure is quite close to the value commonly- 
used in QCD sum rules [33] : 



The above result indicates that our approximation gives a reliable estimative 
of the vacuum energy and that the inclusion of ghosts possibly do not modify 
the value of the vacuum energy [26]. 



6 Conclusion 

We computed the SDE for the gluon propagator in the Landau gauge within 
the Mandelstam approximation where the fermions and ghost fields are ne- 
glected, and where the full gluon vertex is inversely proportional to the gluon 
renormalization function. 

The full triple gluon vertex is also extended to include the possibility of dy- 
namical mass generation, according to a prescription formulated by Cornwall 
many years ago. This prescription does not modify the SDE but is responsible 
for the compatibility of a massive gluon propagator with the Slavnov-Taylor 
identity. 

The renormahzation of the SDE follows a procedure similar to the one pro- 
posed by Cornwall in the renormalization of the SDE in the light cone gauge. 
It is particularly suited for a massive case and leads to a renormalization 
constant of the form = 1 + f{lJ,^, A^). 

We were able to obtain a numerical solution in the full range of momenta that 
we have considered without the need of introducing any asymptotic expression 
for the solution. The propagator is renormalized using the central value of the 
coupling constant at the b quark mass. The solution has been checked to be 
stable within a momentum range of twelve orders of magnitude. 

We verified that the ratio between the dynamical gluon mass and the QCD 
scale {m/AqcD) up to a large extent is independent on the choice of the 
renormalization point, and its value [m/AqcD ~ 2 — 3) is consistent with 
previous estimates for this mass [5]. 

Using a fit to the numerical solution we computed the vacuum energy and 
associated it with the gluon condensate. The value that we obtain is consistent 
with the one usually assumed in QCD sum rules. 




(42) 
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Our calculation is far from being complete and the most obvious extension is 
the introduction of ghosts. An analysis of this case shows that the behavior 
of the gluon propagator is not modified by the inclusion of the ghosts fields. 
However the behavior of the running coupling constant as ^ may be 
different from zero as happens in the present case [26]. 
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